Spreading and shortest paths in systems with sparse long-range connections. 



Cristian F. Moukarzel 
Instituto de Fisica, Universidade Federal Fluminense, CEP 24210-340, Niteroi, RJ, Brazil 



On 



c3 . 



(N 



C 
i 

C/3 



I 

C 

O 
o 



> 

(N 

m 
in 
o 

On 
c3 



c 

o 
o 



X 



Spreading according to simple rules (e.g. of fire or diseases) 
and shortest-path distances are studied on d- dimensional sys- 
tems with a small density p per site of long-range connec- 
tions ( "Small- World" lattices). The volume V(t) covered 
by the spreading quantity on an infinite system is exactly 
calculated in all dimensions. We find that V(t) grows ini- 
tially as F d t d /d for t « t* = {2pT d (d - l)\)~ 1/d and later 
exponentially for t » t* , generalizing a previous result 
in one dimension. Using the properties of V(t), the aver- 
age shortest-path distance £{r) can be calculated as a func- 
tion of Euclidean distance r. It is found that £(r) ~ r for 
r < r c = (2pT d {d - iy.)- 1/d log{2 P r d L d ) and l{r) ~ r c for 
r > r c . The characteristic length r c , which governs the be- 
havior of shortest-path lengths, diverges with system size for 
all p > 0. Therefore the mean separation s ~ p~ x ' d be- 
tween shortcut-ends is not a relevant internal length-scale for 
shortest-path lengths. We notice however that the globally 
averaged shortest-path length t divided by L is a function of 
L/s only. 

PACS numbers: 05.10.-a, 05.40.-a, 05.50.+q, 87.18.Sn 



Regular d-dimensional lattices with a small density 
p per site of long-ranged bonds (or "small-world" net- 
works) jy model the effect of weak unstructured (mean- 
field) interactions in a system where the dominant inter- 
actions have a regular d-dimensional structure, and thus 
may have many applications in physics as well as in other 
sciences p)-p^|. 

In this work we study the spreading (See e.g. p],pT] and 
references therein) of some influence (e.g. a forest fire, or 
an infectious disease) according to the following simple 
law: we assume that, at each time-step, the fire or disease 
propagates from an burnt (or infected) site to all unburnt 
(uninfected) sites connected to it by a link. Long-range 
connections, or shortcuts represent sparks that start new 
fires far away from the original front, or, in the disease- 
spreading case, people who when first infected move to 
a random location amongst the non-infected population. 
For the dynamics of this simple problem, an important 
network property is the set of shortest-path distances 
{iij}, where £^ is defined as the minimum number of 
links one has to traverse between i and j. On isotropic 
c?-dimensional lattices £y is proportional to dfj, the Eu- 
clidean distance between i and j. On regular lattices, 
both the number of sites within an Euclidean distance r 
from i, and the number of sites within r nearest- neighbor 
steps from i behave as r d . 

Disorder can modify this in several ways. On a ran- 
dom fractal, the number of sites contained in a volume 



of Euclidean radius r is r df |1J,|I5[, where df < d is the 
fractal dimension. On the other hand, the number of 
sites visited in at most £ steps is £ d t/ d ™ i ™ |r^ , |l5| , where 
dmin > 1 is the shortest-path dimension, defined such 
that I ~ r dmin . On fractals, shortest-path lengths £{r) 
are thus much larger than Euclidean distances r. 

Consider now a randomly connected network. If we 
have L sites sitting on a regular d-dimensional lattice, 
but connect them at random with an average coordina- 
tion number C (i.e. a total of L d C/2 bonds), the number 
of sites in a volume of radius r is still r d , but we can visit 
~ C k sites in k steps. Thus all L d sites can be visited in 
0(logL d ) steps, and therefore the typical shortest-path 
distance I is of order log L, much shorter than the typical 
Euclidean distance r ~ L. 

"Small- World" networks M are intermediate between 
the regular lattice, where £ ~ L, and the random graph, 
where £ ~ logL. They consist of a regular ci-dimensional 
lattice with N = L d sites, on which pL d additional long- 
range bonds have been connected between randomly cho- 
sen pairs of sites. The key finding of Watts and Stro- 
gatz is that a vanishingly small density p of long-range 
bonds is enough to make shortest-path distances propor- 
tional to logL instead of L. If L d p << 1, the system 
typically contains no shortcuts, and the average shortest- 
path distance £ — l/N 2 J2<ij> [^ij] p scales as L. If on 
the other hand L d p » 1, one finds I - \ogL 0JI0- For 
any fixed density p of long-ranged bonds, a crossover- 
size L*(p) exists above which shortest-path dis- 
tances are only logarithmically increasing with L. This 
crossover size diverges as p~ Y l d when p — » 0. The 
nature of the small- world transition at p = was recently 
discussed controversially j7|,p pl|JI^ ] . It is still a matter 
of debate whether p = can be regarded as being 
equivalent to a critical point. 

In this work we calculate the volume V(t) that is cov- 
ered, on a small- world network, by a spreading quantity 
as a function of time jll] when the spreading law is the 
simple rule above, and derive an exact expression for the 
average shortest-path £(r) on these systems. 

Assume a disease spreads with constant radial velocity 
v = 1 from an original infection site A, as shown in Fig. [I]. 
Let p = 2p be the density of shortcut- ends on the system. 
We work on the continuum for simplicity, so that the 
infected volume will initially grow as a sphere of radius 
t and surface T^f 1 ^ 1 . We call the sphere stemming from 
A "primary sphere" . 

Each time the primary sphere hits a shortcut end, 
which happens with probability pTdt^ 1 per unit time, a 
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new sphere ( "secondary" ) starts to grow from a random 
point in non-infected space (the other end of the short- 
cut). These in turn later give rise to further secondary 
spheres in the same fashion. 




FIG. 1. Consider the spreading of fire or diseases on 
small-world systems. Assume that pL d points are drawn at 
random in d-dimensional space (open dots). The mean dis- 
tance S between neighboring points is proportional to p~ x l d . 
Now connect pairs of points at random. The mean separa- 
tion between "mates" will be of the order of L, the system 
size. Paired points represent long-range bonds (shaded lines), 
across which hre or diseases travel instantaneously. Now a 
disease starts to spread from A. Each time that the resulting 
sphere hits a shortcut-end, a secondary sphere will be born 
somewhere else. These in turn later give rise to other sec- 
ondary spheres. The proliferation of secondary spheres pro- 
duces an exponentially fast growth of the infected volume for 
times t > S. 

Following Newman and Watts [H], we notice that the 
total infected volume is the sum of the primary volume 
T d J Q T d ~ 1 dr plus a contribution V(t — r) for each new 
sphere born at time r. Thus in the continuum the average 
total infected volume satisfies 

V(t) = T d [ r^il + pVit-T^dr, (1) 
Jo 

which can be rewritten in terms of rescaled variables V — 
pV and t = (pT d {d-l)\)^ d t as 

m = jj^j\t-f)^{l + V { f)}df (2) 

It is interesting to notice that V is the total number 
of infected shortcut-ends, while t d /d\ is the total num- 
ber of shortcut-ends infected by the primary sphere. On 
an infinite system, the functional relation (Q) that links 
these two variables has no parameters except for the 



space dimensionality d. On a system of hnite volume 
T d L d /d, an important parameter is the rescaled linear 
size L = (pT d (d - l)\f/ d L, whose d-th. power gives the 
total number N s of shortcut-ends as N s — L d /d\. 
Deriving (|^) d times with respect to t we obtain 

^ d V(t) = l + V(t), (3) 
whose solution is 

Notice that (Q) is a series expansion of (e* — 1) with 
all powers not multiples of d removed. Thus (||) can be 
written as a sum of d exponentials, each with a different 
(i-root of 1 in its argument. In this way, powers which 
are not multiples of d cancel out. 

V(t) = -]Texp{^t}-l (5) 

n=0 

where p d = e l27T ^ d . Some specific examples are 
V{t) = e* - 1 in Id, 

V(i) = coshi - 1 in 2d, 

V(t) = e!+e^!+e^_ 1 m 3d) 

where the one-dimensional solution coincides with that 
previously derived by other methods JO]. 

A general property of (^) is that V grows as i d /d\ 
for t < 1, and later exponentially as e t /d. Thus the 
characteristic timescale |ll| for the spreading process is 

r = (pT d (d-i)ir 1 / d . 

Notice that (0), and thus also (f|), only hold on an 
infinite system. On a finite system with L d >> 1, V 
will saturate after a time t sat that can be estimated by 
equating V ~ e tsat jd ~ L d / dl and therefore 

t sat ~log(L d /(d-l)\), (6) 

which can be rewritten as 

tsat - (pF d (d ~ l)!)- 1 ^ \og( P T d L) (7) 

If on the other hand L d << 1, the spreading stops at 
tsat — L, before reaching the exponential growth regime. 
Thus for a finite system with L d > > 1 one has 

!t d /d\ for t « 1 
e l /d for l«£<t sot ~l g( I3 ^ T ) (8) 
L d /d\ for t>t sa t 
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Assume now that L d » 1. Because of the exponen- 
tially fast spreading process, the fraction of the total vol- 
ume covered by the disease is negligible for t < t sat and 
saturates to one abruptly at t = t sat . Therefore on a large 
system most of the points become infected essentially at 
the same time t sat . 

Now let us see how to calculate the average shortest- 
path distance £(r) as a function of the Euclidean sep- 
aration r between two points. Since we assumed that 
the disease spreads with unit velocity, it is clear that the 
time t at which a point x becomes first infected is exactly 
the shortest-path distance £(A, x) from A to x. By defini- 
tion, no part of the finite system remains uninfected after 
t = t satl so we conclude that no shortest-path distance 
can be larger than t sa t on a finite system. Combining this 
with the fact that £(r) cannot decrease with increasing 
r, we conclude that 



£(r) = t sat 



for r > £, 



(9) 



In order to calculate £(r) for r < t sat , let us write 
V(t) — Vi(t) + V 2 (t), where V\ is the primary volume 
and Vi the volume infected by secondary spheres. 

Let P2(t) be the probability to become infected by the 
secondary infection exactly at time t < t sat . Conse- 
quently I 2 (t) = J Q P2{t)<1t is the probability for a point 
to become infected at time t or earlier. Assuming that 
p 2 (t) is known, it is easy to calculate the average shortest- 
path distance £(r) as a function of Euclidean distance r, 
according to the following. If an individual at x becomes 
infected by a secondary sphere at time r < d E (A, x), its 
shortest-path distance £(A, x) to A is r. Otherwise if x is 
still uninfected at time t = d E (A, x) (which happens with 
probability 1 - I 2 (d E (A, »))), then £(A, x) = d E (A,x), 
since at that time the primary sphere hits x with proba- 
bility one. Therefore the average shortest-path satisfies 



£{r)= / tp 2 (t)dt + r{l-I 2 (r)} 
Jo 



= r - I I 2 {t)dt 
'o 



(10) 
(11) 



The fact that the secondary volume V 2 is randomly 
distributed in space makes this problem relatively simple. 
The probability I 2 {t) for a point to be infected by the 
secondary version of the disease at time t or earlier is 
simply I 2 — V 2 (t) / '(1 / '<fT d L d ) , i.e. the fraction of the 
total volume which is covered by the secondary infection. 
Thus 



£{r) 



r d L d 



V 2 (t)dt 



(12) 



If there are no shortcuts on the system, V2 is zero at all 
times and thus £{r) = r as expected. But it is also clear 
from this expression that £(r) = r when L — > 00, for all 
finite r, i.e. in the thermodynamic limit the shortest-pat 



distance £(r) coincides with the Euclidean distance r for 
all finite r, no matter what p is. 

On a finite system with L d >> 1, V 2 {t)/L d is negligible 
for all t < t sa t a s we have already noticed. Therefore 
£(r) = r if r < t sat . Combining this with (JJ) we have 



£{r) 



r for r < r c 
r c for r > r c 



(pT d (d-l)\)- 1 / d log(pT d L d ) 

(13) 



Detailed knowledge of £(r) for r — r c would only be pos- 
sible if the finite-size effects that we ignored in (||) were 
exactly known, but the interesting remark is that the lack 
of this knowledge has little or no importance for r r c . 

We thus see that on a finite system, a characteristic 
length r c = { P T d {d-l)l)- 1 / d \og(pT d L d ) exists, that gov- 
erns the behavior of average shortest-path distances as 
a function of Euclidean separation. This characteristic 
length diverges when L — > 00, for any p > 0. The typical 
separation s = p~ l l d between shortcut ends which 
is size-independent, is not relevant for £{f). The validity 
of (|l^) has been verified numerically in one dimension 
recently |l3| . 

It is interesting to notice that the rescaled shortest- 
path distance £ — (pT d ) 1 / d £ is a simple function of the 
rescaled Euclidean distance f — (pr ( j) 1 ' d r. 



£(r) 



f for f < 1 
1 for f > 1 



(14) 



Using (|13|) we can now calculate £{p, L) , the (global) 
average shortest-path length d^Q, when L d >> 1. One 
has 



2{p,L) = -t j L £(r)r d - l dr 



1 



L d 
1 /r c \ d 



r d ~ l dr 



d+1 \L 



So that the "order parameter" C — £/ L |9Ul3j] reads 

C = < l ~d^ 



(15) 
(16) 

(17) 



where z = r c /L. 

When p —> faster than L x l d (so that L d « 1), 
formula jl5| ) holds with r c — > L, and thus C — > d/(d + 1) 
as expected. On the other hand if p > one has that 
r c << L when L — > 00, and thus C — > in this limit. 
Therefore C undergoes a discontinuity at p = in the 
L -> 00 limit |Q. 

Notice that r c /L = \og(pT d L d ) / L(pT d (d - l)\f/ d ~ 
log(L / s) / (L / s) , , where s ~ p^ 1 ^ |ll[ is the mean sep- 
aration between shortcut-ends. Thus L can be written 
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as a function of L/s only. Therefore if we measure £ on 
systems with several values of L and p and plot the data 
versus L/s, we would find that they collapse BM^O . 
Because of this behavior some authors have sug- 

gested that the transition at p = is a critical point with 
a size-independent characteristic length £ ~ s ~ ( o~ 1 / d . 
Our results here and in previous work ||[l3| suggest that 
this is not the case. According to our calculation, the 
only characteristic length in regard to shortest-paths is 
r c , and it diverges with system size L. 

We have thus shown that, on a finite system with 
L d p >> 1, two widely separated timescales for spreading 
can be identified. The first one t* = (2pT d (d - 
determines the crossover from normal (i.e. proportional 
to t d ) to exponential spreading. A much larger timescale 
t S at given by (0) determines the saturation of the spread- 
ing process. This second timescale coincides with the 
lcngthscale r c at which the behavior of shortest path 
lengths £(r) saturates, as given by Eq. (|l3|). 

It is clear from our calculation that r c diverges with 
L because the locations of the secondary spheres are un- 
corrected with the location of the primary infection. In 
other words, because on a system of size L, the typical 
separation between both ends of a shortcut scales as L. A 
different situation would certainly arise if shortcuts had 
a length-dependent distribution. For example one can 
connect each site i, with probability p, to a single other 
site j, chosen with probability r/- a , where a is a free pa- 
rameter. For a — > 0, this model is the same as discussed 
here, while for a large one would only have short-range 
connections and thus there would be no short-distance 
regime, even for p = 1. We are presently studying this 
general model fll6[. 
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